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Abstract 

CM , 

We discuss the problem of counting the maximum number of distinct states that an 
isolated physical system can pass through in a given period of time — its maximum 
speed of dynamical evolution. Previous analyses have given bounds in terms of AE, 
the standard deviation of the energy of the system; here we give a strict bound that 
, depends only on E — Eq, the system's average energy minus its ground state energy. 

0^ | We also discuss bounds on information processing rates implied by our bound on 

the speed of dynamical evolution. For example, adding one Joule of energy to a 
given computer can never increase its processing rate by more than about 3 x 10 33 
operations per second. 
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1 Introduction 



In the realm of computation, the first two quantitative questions that one is 
likely to ask about a machine are (i) How much memory does it have? and 
(ii) How fast does it run? In exploring the computational limits of physical 
dynamics, one might try to ask the same questions about an arbitrary physical 
system. 

Question (i) essentially asks, "How many distinct states can my system be put 
into, subject to whatever physical constraints I have?" This is really a very 
old question: the correct counting of physical states is the problem that led to 
the introduction of Planck's constant into physics [11], and is the basis of all of 
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quantum statistical mechanics. This question can be answered by a detailed 
quantum mechanical counting of distinct (mutually orthogonal) states. It can 
also be well approximated in the macroscopic limit [9,21] by simply computing 
the volume of phase space accessible to the system, in units where Planck's 
constant is 1. 

Question (ii) will be the focus of this paper. This question can be asked 
with various levels of sophistication. Here we will discuss a particularly simple 
measure of speed: the maximum number of distinct states that the system can 
pass through, per unit of time. For a classical computer, this would correspond 
to the maximum number of operations per second. 

For a quantum system, the notion of distinct states is well defined: two states 
are distinct if they are orthogonal. The connection between orthogonality and 
rate of information processing has previously been discussed[12,2,7,13,14,4], 
but no universal bound on computation rate was proposed. The minimum time 
needed for a quantum system to pass from one orthogonal state to another has 
also previously been characterized, in terms of the standard deviation of the 
energy A_E[10,16,20,17,5]. This bound places no limit, however, on how fast 
a system with bounded average energy can evolve (since AE can be arbitrar- 
ily large with fixed E). Bounds based directly on the average energy E have 
previously been proposed[6,3], but these bounds apply to the rate of commu- 
nication of bits, rather than to the rate of orthogonal evolution: difficulties 
associated with such bit-related bounds are discussed in [12] . The new bounds 
derived in this paper are also based on average energy, but they apply to rates 
of orthogonal evolution. For an ordinary macroscopic system, these bounds are 
achievable: we show that adding energy increases the maximum rate at which 
such a system can pass through a sequence of mutually orthogonal states by 
a proportionate amount. 

There has recently been much interest in the possibilities of quantum comput- 
ers: computers that can operate on superpositions of computational states[19]. 
Even isolated quantum computers will, in general, pass through sequences of 
(nearly) mutually orthogonal states in the course of their complicated time 
evolutions. At the least, an efficient quantum computation should, for some 
initial states, have a final state that is reasonably distinct from its initial state. 
Thus a bound on the rate of orthogonal evolution is relevant in this case as 
well. 
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2 Maximum rate of dynamics 



In the energy basis, quantum time 
positions of frequency components, 
a maximum energy eigenvalue, the 
should be bounded by 



evolutions are constructed out of super- 
One might expect from this that, given 
frequency with which states can change 



If we take our zero of energy at the ground state of the system,^] and consider 
long evolutions, then this relation is true, as we will discuss below. We will 
also show that, given a fixed average energy E (rather than a fixed maximum), 
there is a similar bound 

IE 

v± < -p (2) 



where again we have taken our zero of energy at the ground state. This equa- 
tion has the following interpretation: in appropriate units, the average energy 
of a macroscopic system is equal to the maximum number of orthogonal states 
that the system can pass through per unit of time. This is the maximum rate 
that can be sustained for a long evolution — the rate at which a system can 
oscillate between two states is twice as great. 



2.1 Orthogonality time 



We begin our analysis by discussing the question of the minimum time needed 
for any state of a given physical system to evolve into an orthogonal state. 

An arbitrary quantum state can be written as a superposition of energy eigen- 
states 

|^ )=E C nl^n) ( 3 ) 



Here and throughout this paper we assume that our system has a discrete 
spectrum, and that the states have been numbered so that the energy eigen- 
values {En} associated with the states {\E n )} are non- decreasing. To simplify 
formulas, we will choose our zero of energy so that Eq = 0. 

2 In all of this discussion, if for some reason no energies below some E m - m are 
allowed, we should take our zero of energy there instead. 
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Let t_i_ be the time it takes for j^o) to evolve into an orthogonal state. We will 
now show that, with a fixed average energy E, it is always true that 

h 

* £ 41' (4) 

This result is somewhat surprising, since earlier results gave a bound only in 
terms of AE 

T - £ ike < 5 > 



This earlier bound would suggest that, given a fixed average energy, one could 
construct a state with a very large AE in order to achieve an arbitrarily short 
t_i_. We will show that this is not the case. 

Let us begin by observing that if j^o) is evolved for a time t it becomes 

\A) =£c n e-^|K> (6) 

n 

Now we let 

oo 

S{t) = falfr) = £ |c n |V^ (7) 

n=0 

We want to find the smallest value of t such that S(t)—0. To do this, we note 
that 



2 / E„t s 



Re(5) = l c «l cos (^ 

71=0 11 

>EN (!--(— + -(—)) 

IK 2 

= l-—t + -1m(S), (8) 
iin 7i 

where we have used the inequality cos a; > 1 — -{x + sinx), valid for x > 0. 
But for any value of t for which S(t) = 0, both Re(S') = and Im(S') = 0, 
and so Eq. (8) becomes 

Thus the earliest that S(t) can possibly equal zero is when t = h/AE, which 
proves Eq. (4). Of course Eq. (8) also gives approximately the same bound on 
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how quickly we can have approximate orthogonality, since if \S(t)\ is small, 
then so are Re(S') and Im(S'). 

This bound is achievable if the spectrum of energies includes the energy 2E 
(and is very nearly achievable if the spectrum includes a value very close to 
this, as we would expect, for example, for any ordinary macroscopic system). 
In this case, we let 

which has average energy E. This evolves in a time t = h/AE into 

w - « (id 

which is orthogonal to l^o)- If we evolve for the same interval again, we will 
be back to j^o): the evolution oscillates between these two orthogonal states. 
For these states, AE = E, and so both of the bounds Eq. (4) and Eq. (5) are 
achieved. 

There are also cases where Eq. (4) gives a much better bound than Eq. (5). 
Consider, for example, the state 

|Vo) = a(|0) + |e)) + b(\ne) + \(n + l)e)) (12) 

which evolves into an orthogonal state in a time r± = h/2e. Given E, as long 
as we choose e < 2E (i.e., r± > h/AE) the average energy of the first pair of 
kets will be less than E. Given e, for large enough n the average energy of the 
second pair of kets will be greater than E. Then we can always find coefficients 
a and b that make the average energy of |^ ) be E and also normalize the state. 
But this state has a AE that depends on our choice of n: in fact AE = Q(y/n). 
With fixed E, AE can be as large as we like. Thus in this case, Eq. (5) is not 
a useful bound and Eq. (4) is optimal. 

2.2 Cycles of orthogonal states 

In the discussion above, we have seen that a quantum system with average 
energy E can be made to oscillate between two orthogonal states with a fre- 
quency of 4E/h. Now we address the question of how fast a quantum system 
can run through a long sequence of mutually orthogonal states. We begin by 
considering the very restricted set of evolutions that pass through an exact 
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cycle of N mutually orthogonal states at a constant rate. In this case it is easy 
to show (see Appendix A) that 



T± > 



N-l h 
N 2E 



(13) 



Thus for very long evolutions that form a closed cycle, the maximum transition 
rate between orthogonal states is only half as great as it is for an oscillation 
between two states. In the next section, we will show that this long-sequence 
asymptotic rate is achievable in principle for any ordinary macroscopic system. 
Here we will first give an example of a system for which an exact cycle of iV 
mutually orthogonal states (cf. [15]) achieves this bound. 

The one-dimensional harmonic oscillator has an exact cycle after some period 
r. Taking our ground-state energy to be zero, all of the energy eigenvalues are 
multiples of e\ = h/r. Let 



If our system passes through N mutually orthogonal states in time r, then 
the average time to pass between consecutive orthogonal states is r± = r/N. 
Noting that eiT±/K = 2n/N, we see that the state obtained from \ip ) after m 
time intervals of length r± is 



N-l 



1 




nei) 



(14) 




(15) 



and so 




n=0 



(16) 



Now we can calculate the relationship between E and t±. 




and so 



r± = 



N-l h 
N 2E 



(18) 
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2.3 Long sequences of orthogonal states 



Now we turn to the question of whether ordinary macroscopic physical systems 
can also run through long sequences of mutually orthogonal states with r± = 
h/2E. We will show by construction that they can. As in the discussion above, 
we will not need to use arbitrarily large eigenvalues to achieve this rate, and 
so our state can be written 

N-l 

|Vo) = E C n\En) (19) 
n=0 



Now we simply let 



IE, 



n+l 

~E 



E n 



(20) 



N 



This definition of c n generalizes our example from the previous section: for 
the special case of E n = ne 1 , Eq. (19) reduces to Eq. (14), which achieves 
t_i_ = h/2E in the macroscopic limit. Notice also that, with this definition of c n , 
states with degenerate energy eigenvalues are not repeated in our superposition 
(they get a coefficient of zero, since the E n 's are numbered in non-decreasing 
order). This definition of c n always gives normalized states, since 

N-l rp rp 

(Wo> = E n % - = 1 (21) 

n=0 E N 

We can calculate the average energy in the state |^ )- This is just 

N-l 771 _ rp 

(MH\^) = E <,c n E n (E n ,\E n ) = E -E n (22) 

n,n' n=0 



For ^> 1 and c n <^ 1, we can approximate this sum by an integral. Letting 

x = n/N and e(x) = E n /E N , we have 



(MH^o) ~ En J e^-dx = ^J j-eHx = ^ (23) 



In Appendix B we estimate the corrections to this approximation, which vanish 
for large N. Thus, with this definition of c n , by giving equal weight to equal 
energy intervals we get an average energy that is half of the maximum energy, 
just as in the limiting case considered in Eq. (17). 
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Now the state obtained from \tp ) after m intervals of length rj_ = h/2E = 
h/E N is 



\^ m ) = Y,^e- 2mm ^\E n ) (24) 



and so 



(VwlVU = E c* n ,cJ^ m ' E «'- mE "\E n ,\E n ) 



n,n' 



E 1 E n+ i — E n 2m^(m'-m) 
e e n k (25) 

n=0 ^ 



As we've seen, this is equal to 1 if w! — m. Let k = m! — m ^ 0, and again 
let x = n/N and e(x) = E n /E N . Then 

\ Jp \. d 2-irike 

(VwlVU « _/ ^e 2mfe£ dx = y ^-^^ = (26) 
o o 



A more careful analysis (see Appendix B) verifies that the corrections to this 
approximate calculation vanish for large N . Thus we can run through a long 
sequence of nearly orthogonal states at the rate u± = 2E/h = E m3iX /h. 



3 Interpretation 

For an isolated macroscopic system s with average energy E^ s \ we have seen 
that we can construct a state that evolves at a rate i/± = 2E^/h. If we 
had many non-interacting macroscopic subsystems, we would have an aver- 
age energy for the combined system of E tot = Y, s E^ ■ Our construction of 
the previous section applies prefectly well to such a composite system, and 
in particular lets us construct a state for this combination of non-interacting 
subsystems that evolves at a rate of is± = 2E tot /h = J2 S 2E^ /h = J^s^i^- 
Thus if we subdivide our total energy between separate subsystems, the max- 
imum number of orthogonal states per unit time for the combined system is 
just the sum of the maximum number for each subsystem taken separately. 
This is analogous to the case in a parallel computer, where the total number 
of operations per second for the whole machine is just the sum of the number 
of operations per second performed by the various pieces. Our result should be 
interpreted in a similar manner: average energy tells us the maximum possible 
rate at which distinct changes can occur in our system. 
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It is interesting to ask how this connection between energy and maximum 
possible number of changes looks in a semi-classical limit. As a simple example, 
let us consider a single-speed lattice gas [8]. This is a classical gas model in 
which identical particles are distributed on a regular lattice. Each particle 
moves from one lattice site to an adjacent site in time ST. At the end of each 
ST interval, all energy is kinetic, and all particles have the same energy SE. 
Thus if the total energy is E, then E /SE is equal to the number of particles, 
and so the maximum number of spots that can change per unit of time is 
equal to 2E/SE: E/SE spots can be vacated, and E/SE new spots occupied. 
Fewer spots will change if some particles move to spots that were previously 
occupied, but we can never have more than 2E/SE changes in time ST. Thus 
if we impose the constraint on this lattice system that SEST > h, our bound 
on the rate at which spots can change becomes z/ changc < IE j SEST < 2E/h. 

It is also interesting to ask how this connection between energy and orthogonal 
evolution looks in different inertial reference frames. Clearly we will see some 
orthogonal states that are purely an artifact of our choice of reference frame. 
For example, an isolated stationary atom in an exact energy eigenstate never 
transitions to an orthogonal state, but if we view this same atom from a moving 
frame we will see a sequence of distinct position states. If we are interested 
in a bound on "useful" dynamics (e.g., on computation), then we shouldn't 
count extra states that arise just from the state of motion of the observer. The 
obvious solution is to define the amount of "useful" dynamics to be the number 
of orthogonal states that the system passes through in its rest frame (center of 
momentum frame). As long as our non-relativistic analysis is valid in the rest 
frame, we can infer that (in that frame) a system with total relativistic energy 
E r cannot pass through more than 2E r t r /h different orthogonal states in time 
t r . Then in any frame, we can compute this bound on "useful" evolution, since 
E r t r = p^x^ . In a frame in which the system starts at the origin at time t — 
and moves in the positive x direction with a momentum of magnitude p, our 
bound is (2/h)(Et — px): from the time component of the bound we subtract 
a space component. Note that we subtract one orthogonal state for each shift 
of a distance h/2p (cf. [5,22]). 



4 Conclusion 

The average energy E (above the ground state) of an isolated physical system 
tells us directly the maximum number of mutually orthogonal states that the 
system can pass through per unit of time. Thus if the system is a computer, 
this quantity is a simple physical measure of its maximum possible processing 
rate. 

Equivalently we can say that Et counts orthogonal states in time. Just as 
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accessible phase-space volume tells us the number of distinct states that a 
macroscopic physical system can be put into, Et "action volume" tells us the 
number of distinct states that a system with energy E can pass through in 
time t. 
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A Minimum orthogonality time for exact cycles 

If we demand that the evolution should exactly cycle after N steps, this condi- 
tion puts a severe restriction on the energy eigenf unctions that can contribute 
to our initial state. If \Ej) contributes to \ipo), and if our cycle length is r, 
then we must have that Ejr/h = 2nkj for some integer kj, which means that 
Ej = kjEi, where S\ = h/r. Thus our initial state must have the form 



For simplicity, we have not included degenerate energy eigenf unctions in our 
superposition — adding these would not affect our conclusions. Thus the as- 
sumption of exact periodicity restricts us to systems with energy eigenvalues 
that are a subset of a one dimensional harmonic oscillator spectrum, in which 
all energies are multiples of h/r. 

As in Section 2.2, we see that the state obtained from |t/> ) after m time 
intervals of length t± = t/N is 



oo 



|-0o> = c n\nei) 



(A.l) 



n=0 



■oo 



c n e 



2-Rinm 
N 



) 



(A.2) 



?1=0 
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and so 



(^|^) = Elcn| 2 e 2 ^ (m '- m) (A.3) 

n=0 



There are actually only A" distinct values of the exponential that appears in 
this sum, since for any integer /, e n (m'-m) _ e 2 7p(m'-m) _ ^y e can co u ec ^ 

i 1 2 1 2 

together all c n 's that multiply each distinct value: let \d n \ = J2^o\ c n+iN\ ■ 
Then 

JV-1 

<Vw|VU = E Mn| 2 e^ (m '- m) = 5 m > m (A.4) 
n=0 



since these states are supposed to be orthogonal. This last equality will ob- 
viously be true if we let |d n | 2 = 1/A" for all ra's. In fact, since there are A" 
different possible values of m! — m, Eq. (A.4) constitutes A" linearly indepen- 
dent equations with A" unknown coefficients, and so this solution is unique. 

i 1 2 i 1 2 

Thus by picking any set of \c n \ 's that add up to make all the \d n \ 's equal to 
1/A", we obtain a state \ipo) that evolves at a constant rate through a sequence 
of A" mutually orthogonal states in time r. 

Now we can calculate the relationship between the average energy and the 
orthogonality time. 

oo N— 1 oo 

E = £i E \cn\ 2 n = ei E E |cn+w| 2 (n + ZiV) (A. 5) 

n=0 n=0 ;=o 



but E/^o |cn+w| 2 (™ + IN) > |c„+«Ar| 2 n = n\d n \ 2 , and since |<i n | 2 = 1/N, 

Since £1 = h/r = h/Nr±, this establishes Eq. (13). Note that we get equality 
in Eq. (13) only for the state given in Eq. (14). 

We can also argue that if we don't use energy eigenvalues larger than -E max = 
ei(N — 1), then there are at most A" eigenstates in our superposition, and so 
we can pass through at most A" different orthogonal states in time r. Thus 

^ T h ( h \ ( N ~ l \ 
r±>- = — -= — — — (A.7) 



N S ± N V^max/ V N 



which yields Eq. (1) for a long exact cycle. 
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B Approximating sums 



Here we estimate the corrections to the average energy and scalar product 
computed in Section 2.3. 

Letting e n = E n j 'En and 5 n = e n+ i — e n , Eq. (22) becomes 



N-l 

(l/j \H\l/j ) =E N J2 ( £ n+l - £ n)tn 
n=0 

E 7V_1 

= E [( £ l+l - £ D - ( £ n+l ~ £ nf 

N-l - 

i - E % 

n=0 

(B.l) 

Now, for any ordinary macroscopic system, from general properties of the 
density of states[18] we know that for large n, E n ~ n c , where c is a positive 
constant much less that one. From this we can show that J2^n — 0(1/ N 2c ). 

Similarly, letting a = 2ix(m! — m) — 2nk, Eq. (25) becomes 




N-l 

n=0 
Y N-l 



ia 



n=0 

N-l oo 



D ia8 n 



-i)-E 



i=2 



(iaS n ) j 



1 JV_i 00 • iaS V 
^_ e ia£n ' 



(B.2) 



n=0 j=2 



Again making the assumption that E n ~ n c for large n, we can show that the 
magnitude of this sum is 0(k/N 2c ). 
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